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0^ ' As a result of climate change, many populations have to modify their 

range to follow the suitable areas - their "climate envelope" - often risk- 
. ing extinction. During this migration process, they may face absolute 

boundaries to dispersal, because of external environmental factors. Con- 
0^ . sequently, not only the position, but also the shape of the climate envelope 

■ can be modified. We use a reaction-diffusion model to analyse the effects 

on population persistence of simultaneous changes in the climate envelope 
position and shape. When the growth term is of logistic type, we show 

■ that extinction and persistence are principally conditioned by the species 
5h I mobility and the speed of climate change, but not by the shape of the 

climate envelope. However, with a growth term taking an Allee effect into 
account, we find a high sensitivity to the variations of the shape of the 
climate envelope. In this case, the species which have a high mobility, 
although they could more easily follow the migration of the climate enve- 
lope, would be at risk of extinction when encountering a local narrowing 
of the boundary geometry. This effect can be attenuated by a progressive 
opening of the available space at the exit of the narrowing, even though 
it transiently leads to a diminished area of the climate envelope. 

Keywords: Alice effect • Biodiversity • Climate change • Climate envelope 
• Conservation • Mobility • Reaction-diffusion • Single species model 
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Introduction 



Over the last century, the global Earth's temperature has increased by about 
0.7°C with, in the past 50 years, an even faster warming trend, of about 0.13°C 
per decade (|IPCQ l 20071) . The consequences of this warming on fauna a nd flora 
are already visible, and well documented (see e.g.. IWalther et al.ll2002l ). Espe- 
cially, poleward and upward shifts of many species ranges have been recorded , 
and are most l ikely linked to the climate change (jParmesan and Yohd l2003t 
|Parmesan|[2006h . 

The vari ations of Earth's climate are in fact highly spatially heterogeneous 
( IPCCII2007D . Thus, depending on the regions they inhabit, species may be more 
or less subject to high changes in climate. Some species can adapt, whereas oth- 
ers, and especially range-restricted sp ecies, risk extin c tion. A recent study based 
on the species-area relationships, by iThomas et all ( 2004 ). predicted between 
15% and 37% of species extinction in the next 50 years due to climate change, 
in sample regions covering 20% of the Earth's terrestrial surface. 

For the next century, climate projections indeed predict a mean increase 
of temperature of 1.8°C for minimum scenarios to 4°C for maximum expected 
scenarios. 

Recently, some authors proposed mathematical models for analysing the 
factors that influence population persistence, when facing a climate change. 
They used the notion of "climate envelope" , corresponding to the environmen- 
tal conditions under which a population can persist. They assumed that the 
conditions defining this envelope did not chan ge with time, while t he envelope 
location moved according to climate change (jThomas et al 1 l2004h . Keeping 
0.9°C/100 km as the poleward temperature gradient, the above-mentioned tem- 
perature increases should imply poleward translations of the climate envelope 
at speeds of 2 — 4.5 km/year. 

A one-dimensional reaction-diffusion model has been proposed bv lBerestvcki et~al 



( 20071 ). Results have been obtained, especially regarding the links between pop- 
ulation persistence and species mobility. Indeed, if the individuals have a too low 
rate of movement, then the population cannot follow the climate envelope and 
thus becomes extinct. On the other hand, if the individuals have a very high rate 
of movement, they disperse outside the climate envelope and thus the population 
also becomes extinct. Other relat ed works, based on one-dimensional reaction- 
diffusi o n models, can be found in Potapov and Lewip (12004 ). Deasi and Nelson 
(|2005l) . |Pachepskv et alj (|2005l) and lLutscher et all (|2006|k 

The model that w e consider here is derived from the cl assical Fisher popula- 
tion dynamics model (|Fisherlll937tlKolmogorov et al.lll937l ). In a two-dimensional 
bounded environment f2, the corresponding equation is: 
du 



dt 



DV 2 u + ug{t, x, u), t e [0, +oo), x G H C 



(1) 



The one-dimensional model considered bv lBerestvcki et all ( 2007 ) is a particular 
case of (p}, with a logistic growth function ug(t,x,u). Indeed, the authors 
assumed the per capita growth rate g to decrease with the population density 
u for all fixed x and t. 
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Other types of growth functions are of interest, especially those taking ac- 
count of Allee effect. Allee effect occurs when the per capita growth rate reaches 
its peak at a strictly positive population density. At low densities, the per 
capita growth rate may then becom e negative ( strong Allee effect). Allee ef- 
fect is known in many species (see lAlleel Il938t iDennid Il989t IVeit and Lewis! 
f996). This results from several processes which can co-occur (IBerec et al. 
20071), such as diminished c hances of finding mates at low densities (|McCarthv 
1997t iRobmet et al.ll2007bft . fitness decrease due to c onsanguinity o r decreased 
visitation rates by pollinators for some plant species ( Groom 19981) . It is com- 
monly accepted that populations sub ject to Allee effect are more extinction 
prone ([Stephens and Sutherland Il999ft . 

In reaction-diffusio n models, Allee effects are generally modelled by equa- 
tions of bistable type ((Fifelll979l iTurchinl H995 IShi and Shivaiill2006ft . Math- 
ematical analyses involving these equations have demonstrated important ef- 
fects of the domain's geometry, especially while studying travelling wave solu- 
tions. These solutions generally describe the invasion of a constant state, for 
instance where no individuals are present, by another constant state (typically 
the carrying capacity), at a c o nstant speed, and with a con stant profile (see 
Aronson and Weinberger 1978 ). Berestvcki and Flame] ( 2006f ) have proved that 
in an infinite environment with hard obstacles, but otherwise homogeneous, 
travelling waves solutions of bistable equations may exist or not, depending 
on the shape of the obstacles. In the s ame i dea, in an infinite homogeneous 
square cylinder, IChapuisat and Grenier (l2005ft proved that travelling wave so- 
lutions may not exist if the cy linder's diameter is suddenly increased somewhere. 
See al so iMatano et al.l ( 2006) for another related work. In parallel, Keitt et al 



( 20011 ). while studying invasion dynamics in spatially discrete environments, 
showed that Allee effect can cause an invasion to fail, and can therefore be a key- 
factor that determines the limits of species ranges. In on e-dimensional models , 



Allee effect has also been shown to slow-down invasions (IHurford et al.l 2006). 



or ev en to stop or rev erse invasions in p resence of predato rs ( Owen and Lewis! 
2001ft , or pathogens (|Hilker et al.ll2005l ). In a study by iTobin et alj (|2007ft . 



empirical evidences have also been given regarding the fact that geographical 
regions with higher Allee thresholds are associated with lower speeds of invasion. 

The aim of this work is to study, for the simple two-dimensional reaction- 
diffusion model ([I]), how the population size variation during a shift of the 
climate envelope depends on the shift speed, the geometry of the environmental 
boundary, and on the population mobility and growth characteristics. 

The first section is dedicated to the precise mathematical formulation of 
the model. Two types of growth functions are considered, logistic-like, or tak- 
ing account of an Allee effect, with in both cases a dependence with respect 
to the climate envelope's position. We define three domain types, correspond- 
ing to three kinds of geometry of the environmental boundaries. Domain 1 is a 
straight rectangle, and domain 2 is the union of two rectangles of same width by 
a narrow corridor. The comparison between domains 1 and 2 enables to analyse 
the effects of a local n arrowing of the habitat on pop ulati on persistence. As sug- 
gested by the work of Chapuisat and Grenier ( 2005ft and Berestvcki and Hamell 
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(2006), in the case with an Allee effect, extinction phenomena may not be caused 
directly by the reduction of climate envelope due to the narrowing of the domain 
in the corridor, but by its too sudden increase at the exit of the corridor. This 
is why we introduced domains of type 3, which correspond to the union of two 
rectangles of same width by a narrow corridor, which gradually opens over a 
trapezoidal region of length h (the case h = corresponds to domain 2) . Results 
of numerical computations of the population size over 30 years are presented 
and analysed, under different hypotheses on the growth rate, mobility, speed of 
climate change, and domain's shape. These results are further discussed in the 
last section of this paper. 



Formulation of the model 

The population dynamics is modelled by the following reaction-diffusion equa- 
tion: 

du o „ „ 

— = DV 2 u + ug(t,x,u), t e [0, +oo), x e Q C M. 

Here, u = u(t, x) corresponds to the population density at time t and position 
x = {x\,X2)- The number D > measures the species mobility, and V 2 stands 

2 ® 2u & 2u 

for the spatial dispersion operator V u = ——^ + t—^t. The set Q is a bounded 

ox\ OX2 

subdomain of R 2 . We assume reflecting boundary conditions (also called no-flux 
or Neumann boundary conditions): 

— (t,x) = 0, for x e on, 

on 

where <9f2 is the domain's boundary and n = n{x) corresponds to the outward 
normal to this boundary. Thus, the boundary of the domain J7, or equivalently 
the environmental boundary, constitutes an absolute barrier that the individuals 
cannot cross. 



Growth functions 

The function g corresponds to the per capita growth rate of the considered 
species. In our model, it can be of two main types, g\ or g a . The first case 
corresponds to a logistic-like growth rate, depending on the position with respect 
to the climate envelope: 

f gi (t,x,u)=r+(l-%), ifzeC(t), 

\ gi (t,x,u)=r+(l-%)-r-, ifx^C(t). [Z) 

The set C(t) C ^ corresponds to the climate envelope at time t, the real number 
r + > is the intrinsic growth rate of the species inside the climate envelope, 
K corresponds to the carrying capacity inside the climate envelope, and r~ 
corresponds to the drop in intrinsic growth rate, outside the climate envelope. 

In the second strong Allee effect is modelled. For fixed values of 

t and x, the function g a (t, x,u) does not attain its maximum at u = 0, and 
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Figure 1: (a): profile of the growth rate u i— > ugi(t, x,u) inside the climate 
envelope (x € C(t), solid line), and outside the climate envelope (x ^ C(t) dashed 
line) in the logistic case, (b): profile of the growth rate ug a (t,x,u) inside the 
climate envelope (solid line), and outside the climate envelope (dashed line), in 
the case with strong Allee effect. We have set r+ = 1, r~ = 2, K = 10, and 
p = 0.25. 



furthermore g a (t, x,0) < for all t > 0, x S ft. The typical form of per capita 
growth term taking account of an Allee effect is 

/ u\ ( u — c 

where r is a growth term, k is the environment's carrying capacity a nd c is the 
Allee threshold (see e.g.. iLewis and Kareivalll993t iKeitt et al.l 120011 ). For the 
comparison with the logistic case @ to stand, we choose r and k such that 
g a (t, x, K) = gi(t, x, K) — for all x € C(t) and 

max g a (t,x,u)— max gi(t, x,u) = r + . 

ue(0,K),x£C(t) u£(0,K),x£C(t) 

Moreover, we impose an "Allee threshold" c equal to pK, for p E (0, 1), which 
means that g a (t, x,u) < as soon as u < pK, and for all x G Q. Such conditions 

yields k = K and r = -rr. As in the logistic case, we assume that g a drops 

(i-p) 2 

by r outside the climate envelope. Finally, we obtain, 

f ga (t,x,u)=r+ JT ±y(l-%)(j t -p) i£x€C(t), 



g a (t,X,u)=r+j£y(l-%)(%-p) 



\ix<£C(t). 



(3) 



The profiles of the growth functions ugi and ug a are depicted in Figure [T] 
for x inside and outside the climate envelope. 

Throughout this paper, we make the hypothesis r + < r~ . 

Remark 1: In both cases, the environment's carrying capacity is equal to K in- 
side the climate envelope. However, outside the envelope, the carrying capacity 
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in the logistic case would be K yL —j for r + > r . However, it is not de- 
fined when r + < r~ . Indeed, in this last situation, the equation ugi(t,x,u) = 
has no positive solution. Similarly, in the case with Allee effect, the equation 
ug a {t,x,u) = has no positive solution outside the climate envelope as soon 
as r + < r~ , since max q a (t,x,u) = max qi(t,x,u) = r + . This 

uE(0,K),xEC(t) u£(0, K),x&C{t) 

means that the environment is not suitable for persistence outside the climate 
envelope. 

As explained in the introduction, we assume that the climate envelope C(t) C 
moves poleward, according to the climate change. Assuming we are in the 
Northern hemisphere, we consider this move to be of constant speed v > 0, in 
the second variable direction (to the "North"). Furthermore, we assume C(t) to 
be of constant thickness L > 0: 

C(t) := {x = {x\,X2) € ^, such that xi £ [vt, L + vt)}. 

Remark 2: In this framework, the spreading speed of the population, if it 
survives, is constrained by the speed at which the climate envelope moves. 
Let us recall that, in an homogeneous one-dimensional environment, with a 
logis tic growth rate ur + (\ — u/K), the spreading speed is c\ = 2 V r + D (see 
e.g., Aronson and Weinbergerlll978| ). In such an homogeneous environment, 
but with the growth rate, 4wr+/(l — p) 2 (1 — u/K) (u/K — p), taking account 
of a n Allee effect, the spreadi ng speed becomes c* = 2-\A'+Z?v / 2(l/2 — p)/(l — p) 
(see Lewis and Kareiva 19931 ). We observe that c* < c*, and that c* decreases 



with p. Thus, the spreading speed is slowed down by the Allee effect. The case 
p = corresponds to weak Allee effect: the per capita growth rate exhibits a 
maximum not at u = 0, but remains positive at low densities. In that case 
c*[p = 0] = c*/\/2. The cases p > 1/2 corresponding to very strong Allee 
effects, lead to population extinction fsee lLewis and Kareiva 1993 ). 



Geometry of the environmental boundary 

We work in piecewise C 1 bounded domains fl of M 2 of three different types: 

1 A straight rectangle. 

2 The union of two rectangles of same width by a narrow passage u>. 

3 The union of two rectangles of same width by a narrow passage u> followed 
by a trapezoidal region of height h. 

These three domain types, and their position relatively to the origin, are 
depicted in Figure 

The domains are assumed to have a width of 20 km (on the larger parts, for 
domains 2 and 3). In domains 2 and 3, the "Southern" rectangles are assumed 
to have a length of 40 km. The width of the passage u> is assumed to be 2 km, 
and its length 4 km. The "Northern" environmental boundaries are situated 
far enough (400 km away from the "Southern" boundaries), so that they have 
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no influence on the population. The climate envelope is assumed to have a 
latitudinal range L — 30 km. 

Initial condition 

We assume that the population was at equilibrium, and that the climate enve- 
lope was stationary, before the considered period starting at a time t — 0. This 
means that the initial condition w(0, a;) is a positive solution of the stationary 
equation: 

- D\7 2 p = p g(0,x,p), xe Slcl 2 , (4) 

with Sj^(x) — over dft. We place ourselves under the appropriate conditions 
for existence of such positive stationary solutions, in the logistic and Allee effect 
cases (see the Appendix for more details). 



Results 

Using a second order finite elements method, we computed the solution u(t, x) of 
the model (Q} with the growth functions and initial conditions discussed above. 
We focus here on the population size, 

P(t) := / u(t,x)dx, 



Jn 

which was computed over 30 years in various situations. 

Unless otherwise mentioned, we set r + = 1 year -1 for the intrinsic growth 
rate coefficient inside the climate envelope, and we assumed that the per capita 
growth rate was decreased by 2 outside the climate envelope: r - = 2 year -1 . 
We assumed that K — 10 individuals/km 2 . In our computations, the diffusion 
coefficient D varied between 1 km 2 /year, corresponding to populations with 
low mobility, e.g., some insects spe cies, and 50 km 2 /year, corresp onding to 
populations with high mobility (see IShigesada and Kawasaki 1997 . for some 



observed values of D, for different species). We first chose p = 0.25, so that 
the Allee threshold pK is 2.5 individuals/km 2 . The speeds used for the climate 
envelope shift varied between v = 1 km/year and v — 6 km/year. 



Population size over time 

In order to get a first insight into the general behaviour of the population size 
P(t) for t G [0, 30], we first computed it with the fixed values v = 2.5 km/year, 
and D = 10 km 2 /year. The results are presented in Figure[3l 

First notice that, as expected by our choice of growth functions, the initial 
population sizes, -P(O), are almost the same with both growth functions @ and 
([3]): 5800 and 6000 respectively in the logistic case and Allee effect case. A 
slight transient increase or decrease of P(t), probably due to the initial effect of 
the Southern boundary, vanishes after 5 years. 

In domain 1, for both growth functions ([2]) and ([3]), after this short period 
of 5 years, the population sizes remain stable around their initial values. 
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Figure 2: Three types of domains SI. The shaded area corresponds to the climate 
envelope C(0) at t = 0; it is represented here with a length L = 30km. The 
narrow passage w of domains 2 and 3 is separated by thin dotted lines, and is 
followed by a trapezoidal region of height h in domain 3. 
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Figure 3: Computation of the populations size P(t) in function of time t. (a) 
Case of the logistic growth rate ugi(t,x,u), in domain 1 (solid line), and in 
domain 2 (dashed line), (b) Case of the growth rate ug a (t,x,u) with Allee 
effect, in domains 1 (solid line), 2 (dashed line) and 3 (dotted line). For these 
computation, we fixed D = 10, v = 2.5, r + = 1, r~ = 2, K — 10, p — 0.25 and 



On the contrary, the populations react very differently in the domain 2, 
depending on the type of their growth rate. In the logistic case (Fig. the 
population size recovers its initial value, after a transient decrease that lasts as 
long as the corridor is included in the climate envelope. In the case with Allee 
effect (Fig. [3Jd), the population size declines to after 20 years. 

The behaviour of the population with Allee effect in domain 3, with h = 
25 km (Fig. \3jh>), leads, as in domain 2 with logistic growth, to the recovery of 
the initial level of population. 

As a preliminary conclusion of this first step, it comes that the fate of the 
population will be driven by the interaction between environmental parameters 
(v and the type of domain) and biological parameters (D, r + , r~ , K, p, and 
the type of growth function). 

Intertwined effects of mobility and environmental parame- 
ters 

We computed the "final" population size P(30), for the range of parameters 
(v, D) e (1, 6) x (1, 50). We define the "extinction region" as the portion of the 
parameter space leading to P(30) < 1. 



h = 25. 
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(a) domain 1 



(b) domain 2 



Figure 4: Logistic growth rate ugi(t, x,u). Contour lines of P(30) in function 
of the speed of climate change v and of the diffusion coefficient D, in domains 
1 (a) and 2 (b). The hatched areas correspond to the extinction regions. We 
fixed r+ = 1, r~ = 2 and K = 10. 



Logistic case 

The extinction region is very similar in domains 1 and 2 (Fig. [4]), and corre- 
sponds to high values of v, and low D values. In those cases, the population 
mobility (D) is not sufficient to follow the climate envelope (moving at speed v). 
In domain 1, for each value of D, P(30) is decreasing with respect to v, while in 
domain 2, P(30) is decreasing in v for v > 1.5 (whenever v < 1.5, the climate 
envelope has not totally cro ssed the corridor ui at t = 30). As already proved 
in the one-dimensional case iBerestvcki et al. I (|2007h . and for the same reasons 



recalled in the introduction section, we here observe that population sizes P(30) 
are not monotonically linked with the parameter D. 



Case with Allee effect 

In domain 1 (Fig. [5^,), P(30) is decreasing with respect to v and increasing 
with D. The extinction region again corresponds to high values of v, and low 
D values. However, the extinction region is wider than in the logistic case. In 
domain 2 (Fig. [SJd), the population does not survive as soon as v > 1.8. The 
extinction region is therefore very wide compared to the one of domain 1. 

Contrarily to the logistic case, a marked difference between the fate of pop- 
ulations in domains 1 and 2 does exist. This reflects the sensitivity to environ- 
mental boundary geometry, of populations subject to an Allee effect. In domain 
3, with h = 25, where the opening at the exit of the corridor is more progressive 
than in domain 2, we observe that the extinction region is reduced (Fig. [5b) . It 
remains however larger than in domain 1, especially for large values of D. 
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Figure 5: Case with Alice effect, (a), (b) and (c): contour lines of P(30) in 
function of v and D in domains 1, 2 and 3, respectively. The hatched areas 
correspond to the extinction regions. For these computations, we fixed r + = 1. 
r~ = 2, K = 10, p = 0.25 and h = 25. 
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Figure 6: Case with Allee effect. Contour lines of P(30) in domain 3, in function 
of D and of the height h of the trapezoidal region. The hatched area corresponds 
to the extinction region. The other parameters values were r + = 1, r~ = 2, 
K = 10, p = 0.25 and v = 2.5. 

Allee effect and boundary geometry 

Effect of gradual boundary change at the exit of the corridor 

We have seen, in the case of population subject to an Allee effect, that do- 
main 3, with h — 25, gives population persistence more chances than domain 2 
(which corresponds to domain 3, with h = 0). Let us fix v — 2.5, and see how 
persistence depends on h and D. The results of the computation of P(30), for 
(h,D) G [0,30] x (1,50) are presented in Figured 

For each value of D, P(30) increases with h; thus higher h values give popu- 
lation persistence more chances. Furthermore, for each D > 3, we observe that 
we can define a real number h* such that the population goes to extinction if 
h < h* , and the population survives if h > h* . Remarkably, h* defines two 
types of geometry, characterised by the progressiveness of the opening at the 
exit of the corridor u>, which lead to extinction or survival, respectively. For 
most values of D {D > 5), h* linearly increases with D. For small values of D 
(D < 3), extinction occurs independently of the value h, which simply means 
that, as in domain 1, the population cannot follow the climate envelope shift. 

These results indicate that the increased extinction risk in domain 2, com- 
pared to domain 1, may not be caused by the shrinking of the climate envelope 
and its fragmentation in two parts, but by the lack of progressiveness of the 
opening at the exit of the corridor. 
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Figure 7: Case with Allee effect. Solid lines: Number of individuals Pi(t) in the 
region Hi — il n {y £ (44,48)}, corresponding to the four kilometres following 
the corridor, in domains 2 and 3. Dashed lines: population density u(t, c) at the 
"centre" x c = (10; 46) of Qi. For this computation, we fixed D = 10, r + = 1, 
r~ = 2, K = 10, p = 0.25 and v = 2.5. 

Population density at the exit of the corridor governs the fate of the 
population 

Let us focus on what happens at the exit of the corridor to. We considered the 
region f^i = £1 n {y S (44,48)}, corresponding to the four kilometres following 



in this region, and the population density u(t,c) at the "centre" x c of tlx, of 
coordinates (10; 46). We compared domains 2 and 3, with the fixed parameters 
values v — 2.5, D — 10, p — 0.25, and h — 25 (which is larger than h* = 10 in 
this case, see Fig. [6]). The results are presented in Figure [71 

We observe that, at the beginning of the invasion of Sli, for t 6 [4, 7], the 
populations P\(t) are almost identical in the two domains. Thus, the boundary 
geometry after the corridor has no effect on the number of individuals which 
first invade the region On the contrary, during the same period of time, the 
population density at x c , which is situated at the exit of the corridor, remains 
lower in domain 2, and never reaches the Allee threshold pK = 2.5. Thus, 
in domain 2, the dispersion of the few invaders into an open wide space, after 
the corridor, leads to a low population density, and therefore to a negative 
value of the per capita growth rate g a (t,x,u) (it would be false for a logistic 
growth function). The higher density in domain 3, for h > h* , results from the 
progressiveness of the domain width offered to population development. Even if 
the individuals are scarce in fii, the population density at x c exceeds the Allee 
threshold at t — 7. That makes the population increase, leading to survival 
(survival can indeed be observed in Fig. [SJ;). Thus, the fate of the population 
seems to depend on whether or not the population density at the exit of the 



the corridor. We computed 




u(t,x)dx of individuals 
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Figure 8: Case with Allee effect. Minimal value h* of h required for persistence, 
in terms of the parameter p. The other parameters values were D = 10, r + = 1, 
r~ = 2, K = 10 and v = 2.5. 

corridor reaches the Allee threshold. 

This leads us to analyse the effect of the Allee threshold. For the same values 
v = 2.5 and D — 10, the function p ^ h* is presented in Fig. [H For small 
values of p {p < 0.16), we observe that h* — 0, i.e., the population even survives 
in domain 2. Then, h* increases with p. For higher values of p, the population 
never persists, at least when h < 30 (see also the end of remark 2). Finally, the 
stronger the Allee effect, the more progressive must be the opening at the exit 
of the corridor lo. 

On the contrary, the degree of hostility outside C(t), quantified by r~, does 
not seem to affect much h* . Indeed, computing h*, for v — 2.5, D = 10 and 
p = 0.25, we always obtained h* = 10, for r~ E (1,2) (remember that, in this 
work, it is assumed that r~ > r + , see remark 1). 

Discussion 

Using a two-dimensional reaction-diffusion model, we have studied the fate of 
populations with different mobility and growth characteristics, facing environ- 
mental changes: a shift of their climate envelope, the shape of which can be 
modified by environmental boundary geometry. 

The growth functions we considered were of two types, logistic, with a per 
capita growth rate modelled by @ or taking account of a strong Alice effect 
([3]). Boundary geometry diversity has been summarised into three schematic 
domain shapes (Fig. [2]). 

In the logistic case, the other biological parameters being fixed, the response 
(survival or extinction) to a climate envelope shift at speed v is simply deter- 
mined by the mobility of the population, which is measured by the diffusion 
parameter D. The minimum required mobility for survival increases with the 
speed v. This response is independent of the environmental boundary geometry. 
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This clearly appeared comparing domains 1 and 2. In fact, we conjecture that 
it is true for any local perturbation of the boundary (it could be proved using 
the methods of Berestycki and Rossi, in press). 

In the case with Allee effect, we observed a more complex pattern of interac- 
tions among biological and environmental parameters. In the straight domain 
1, similarly to the logistic case, the higher the climate envelope shift speed, the 
higher the required mobility for survival. Additionally, this required mobility, 
at any v, is always higher than in the logistic case. 

Facing a local narrowing of the environment, as in domain 2, contrarily to 
the logistic case, the chances of survival dramatically decrease in the presence 
of an Alice effect. In that case, high mobility is not sufficient to face climate 
envelope shifting at high speeds. Indeed, the population density drop at the 
exit of the corridor exhibits the sensitivity of this type of populations to low 
densities. This drop can be attenuated by a progressive opening of the available 
space at the exit of the corridor, as in domain 3. Even though this geometry 
transiently leads to a diminished area of the climate envelope, it finally results 
in higher chances of survival. The opening has to be all the more progressive 
since the Allee threshold and the species mobility are high. 

Thus, a population subject to an Allee effect should have a mobility which 
is neither too low, in order to be able to follow the climate envelope, nor too 
high, in order to overcome various changes of the shape of the climate enve- 
lope. For this reason, species having an increasing dependence of their mobility 
with respect to their d ensity should be more robust to climate change (see 
Okub o and Levinll2002l for a survey of reaction-diffusion models with density- 



dependent dispersal terms). On the other hand, under the environmental changes 
considered in this paper, populations with negatively density-d ependent disper- 
sal, such as in the model proposed bv lKing and McCabel (|2003l ) for the dispersal 
of early Palaeoindian people in North America, should have high probabilities 
of extinction, if an Allee effect was assumed. 

The reflective boundary conditions that we used throughout this paper mean 
that individuals encountering the environmental boundary are reflected inside 
the domain. These boundary conditions can be encou ntered in many real- world 
situa t ions, corresponding to c liffs, rivers or coasts (phigesada and Kawasaki 
1997 ; Jaeger and Fahrig] [2004) . Other boundary conditions could have been 
considered. With Robin boundary conditions for instance, a part of the indi- 
viduals c rosses the boundary. These boundary conditions write du/dn + eu = 
(see e.g., ICantrell and Cosner 2003L for some details on this type of boundary 
conditions). Preliminary numerical computations have shown that our results 
still hold in such a case, at least for small positive values of e, corresponding to 
few individuals cros sing the environm ental boundary. 

It was stated in iParmesanl (|2006r ). on the basis of empirical studies, that 
range-restricted species, like polar and mountain-top species, were at high risk 
of extinction induced by warming. Our study suggests that some other species 
may fail to expand poleward and that their capacity to expand is linked to the ge- 
ometry of the geographical limits. Apart from laboratory tests under controlled 
conditions, the diversity of arrangements of the Alpine valleys could allow us 
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to see whether the results of this paper can be observed in natural conditions. 
Indeed, insects such as pine processionary moth, t he present range exp ansion 



of which is undoubtedly related to climate change (jRobinet et al.ll2007a[) . may 
allow us to test statistically the relationships between insect progression and 
the geometry of the Alpine corridors. 

In this paper, a southern retraction of the climate envelope was assumed. 
This can be directly linked, for some species, to the fact that they are sensitive 
to high temperatures. This can also be caused indirectly by competition with 
other species, the range of which is shifting to the North. For some species, 
however, this retraction does not occur, leading to an expansion of the species 
range. This could be easily integrated in our model, by setting 

C{t) := {x = (xi,x 2 ) G Cl, such that x 2 € [0,L + vt]}. 

In this case, comparable results should be obtained, with stagnation instead of 
extinction. 



Appendix: initial conditions 

When the growth rate is of logistic type, with g — gi satisfying |(5J), a sufficient 
condition for the existence of a positive solution p i of (|H) can b e derived by 



finding an appropriate sub-solution of ((4} (see e.g., Amannl[~L976l ): indeed, the 
constant K is readily a super-solution. With our choices of f2, C(0) and boundary 
conditions, for a > small enough, the function 

(7T Try \ 
- + —j for y < L, and <f>(x, y) = for y > L, 

n 2 

is such a sub-solution, as soon as -jfa- < r + . Under this condition, the function 
Pi is then the unique pos itive and bounded solution of ^) (it can be proved as 
in Berestvcki at al. 2005f) . In the logistic case, we assume that u(0, x) = p;. 



In the case ([3]) with Allee effect, the condition for the existence of a solution 
of ([4]), with g — g a is more complex, and multiple solutions may exist. However, 
the existence of a positive and bounded solution of (J4J) is still granted when 
C(0) contains a sufficiently large ball Br. Indeed, for R large enough, there 
exists a positive solution i\) of ([4]) on Br, with Dirichlet boundary conditions 
(jBerestvcki and Lionslll980h . The function ip, extended by to Q, is then a 



sub-solution of ([U; the constant K is again a super-solution; this implies the 
existence of a solution p a of (|4]) on fl We can assume in this case that u(0, x) = 

Pa- 

For our computations, pi was obtained by numerically solving (0} thanks to 
a second order finite elements method with triangular mesh elements. In the 
case with Allee effect, p a was computed as the limit of the solution u>(T, x), as 
T — > +oo of the initial-value problem 
dw 

— - DV 2 w = w g(0, x, w), t>0, x e fl, 

with w(0,x) = 2K and §^(t, x) — over dQ. The solution w{T,x) was also 
obtained thanks to a finite elements method. The convergence of w(T, x) to a 
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positive equilibrium, as T —* +00, ensures that the condition for the existence 
of p a is fulfilled. 
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